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Abstract 


Environments play an important role in galaxy formation and evolution, particularly in regulating the content of neutral 
gas. However, current H I surveys have limitations in their depth, which prevents them from adequately studying low H I 
content galaxies in high-density regions. In this study, we address this issue by employing the Five-hundred-meter 
Aperture Spherical radio Telescope with extensive integration times to complement the relatively shallow Arecibo 
Legacy Fast Arecibo L-band Feed Array H I survey. This approach allows us to explore the gas content of dwarf galaxies 
across various environments. We observe a positive relationship between H I mass and stellar mass in dwarf galaxies, 
with a well-defined upper boundary for H I mass that holds true in both observations and simulations. Furthermore, we 
find a decrease in the H I-to-stellar mass ratio (My ;/M,,) as the density of the environment increases, irrespective of 
whether it is determined by the proximity to the nearest group or the projected number density. Comparing our 
observations to simulations, we note a steeper slope in the relationship, indicating a gradual gas-stripping process in the 
observational data. Additionally, we find that the scaling relation between the My ;/M,, and optical properties can be 


improved by incorporating galaxy environments. 
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1. Introduction 


The baryons in the Universe are primarily composed of 
hydrogen and helium, with hydrogen making up around 75% of 
the total baryonic content. In the paradigm of the standard galaxy 
formation model (e.g., White & Rees 1978; Springel et al. 2005; 
Guo et al. 2011) galaxies form when the gas cools, condenses, and 
forms stars at the centers of their dark matter halos. When dark 
halos merge, their galaxies are accreted into more massive 
systems and orbit as satellite galaxies in groups and clusters. 

HI gas, serving as the essential fuel for star formation, is 
indispensable for our understanding of galaxy formation and 
evolution. Different from the stellar component, H I gas could 
be more susceptible to the environments as it can extend to 
larger distances compared to optical radii (Haynes et al. 1984). 
A significant portion of galaxies in the local Universe reside in 
groups and clusters (Eke et al. 2006; Robotham et al. 2011), 
where they suffer from various environmental influences such 
as tidal interactions, harassment, ram pressure stripping, and 
evaporation (e.g., Gunn & Gott 1972; Moore et al. 1996; Dénes 
et al. 2016; Cortese et al. 2021; Rhee et al. 2023). 

Over the last decades, numerous ground-based surveys: such as 
2dF Galaxy Redshift Survey (Colless et al. 2001), the Sloan 
Digital Sky Survey (SDSS; York et al. 2000), the United Kingdom 
Infrared Telescope Infrared Deep Sky Survey (Lawrence et al. 
2007), the Galaxy and Mass Assembly (Driver et al. 2011), the 


Two Micron All Sky Survey (MASS; Skrutskie et al. 2006), 
Redshift Survey (2MRS; Huchra et al. 2012), the Dark Energy 
Spectroscopic Instrument (DESI) surveys (Dey et al. 2019; DESI 
Collaboration et al. 2023a, 2023b) and space-based surveys: such 
as the Galaxy Evolution Explorer (Martin et al. 2005), the Wide- 
field Infrared Survey Explorer (Wright et al. 2010), the Cosmic 
Evolution Survey (Scoville et al. 2007), the Cosmic Assembly 
Near-infrared Deep Extragalactic Legacy Survey (Grogin et al. 
2011), and the ongoing James Webb Space Telescope surveys 
have obtained vast amount of photometric and spectroscopic data 
for millions of galaxies in the local Universe and at high redshifts. 
These successful surveys have provided valuable insights into the 
evolution of stellar masses and star formation rates in galaxies (see 
the review paper and the reference in Madau & Dickinson 2014). 
However, to fully understand the process of galaxy formation and 
evolution, it is crucial to understand the mechanisms of gas 
accretion onto galaxies and the efficiency with which that gas is 
converted into stars (e.g., Kennicutt 1998; Kereš et al. 2005; Salim 
et al. 2007; Leroy et al. 2008). 

An economical approach to obtain the HI is via galaxies’ 
optical properties. Kannappan (2004) found that the H I-to- 
stellar mass ratio, My 1/M,,, exhibits a strong correlation with 
optical-optical (e.g.; 4 —r) and optical-NIR (e.g., u — K) 
colors, allowing for its estimation with a typical scatter of 
approximately 0.4 dex. Subsequently, several studies have 
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aimed to enhance the accuracy of photometric estimators for 
My 1/M,,, by including stellar surface mass density, and optical 
(Zhang et al. 2009) or near-ultraviolet-optical (Catinella et al. 
2010; Zhang et al. 2021) colors. These estimators typically 
yield a scatter of around 0.3 dex in log My 1/M., providing 
improved H I-to-stellar mass ratio scaling relations. 

Direct measurement of HI relies on radio observations, 
which recently have also collected unprecedented data. For 
example, the blind HI survey, the 10096 complete Arecibo 
Legacy Fast Arecibo L-band Feed Array (ALFALFA; Haynes 
et al. 2018, o.100) has detected approximately 30,000 
extragalactic HI line sources within a redshift of 0.06. 
However, due to the limited efficiency with an average 
integration time of 48s for each source, this catalog is rather 
shallow, with a flux limit of 0.18 Jy for a typical velocity width 
of 50 kms !, corresponds to 10°”*M., in HI mass at z ~ 
0.0265. The Five-hundred-meter Aperture Spherical radio 
Telescope (FAST; Nan 2006; Jiang et al. 2019), with a 
300 m effective diameter and the ability to move within around 
50°, has the potential to significant expand nearby HI 
surveys. Similar to the ALFALFA blind survey, using the 
FAST Commensal Radio Astronomy FasT Survey (Zhang et al. 
2021) aims to detect nearly 4.8 x 10? galaxies up to a redshift 
of 0.1, covering a sky area of 20,000 deg’. In this work, we use 
FAST to select 14 galaxies randomly from high-density and 
low-density environments and increase the integration time to 
around 19.2 minutes. This enables us to include gas-poor 
galaxies in various densities and shed light on how galaxies 
evolve in different environments. 

In this study, we use a volume-limited sub-sample of the 
largest blind HI survey, the 40% complete ALFALFA 
(a.40; Haynes et al. 2011) catalog and gas-poor galaxies from 
FAST observation to examine the H I-to-stellar mass ratio in 
different environments, and the scaling relations of H I-to- 
stellar mass ratio with optical properties. In Section 2, we 
briefly describe the sample selection criteria, the observation 
details, the data reduction, and the methods used to extract 
relevant physical properties. Our main results are presented and 
discussed in Section 3. We summarize our results in Section 4. 


2. Data and Method 


Our aim is to investigate the H I-to-stellar mass ratio in 
different environments. H I fluxes are obtained using FAST and 
a.40, while stellar mass is obtained using the data from the 
SDSS survey. 


2.1. ALFLAFA Dwarf Galaxies 


We select galaxies from the matched sample between the 
40% complete ALFALFA (a.k.a. a.40; Haynes et al. 2011) 
survey and the Seventh Data Releases of SDSS (SDSS DR7; 
Abazajian et al. 2009). 
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The SDSS uses a 2.5 m telescope located at Apache Point 
Observatory near Sacramento Peak in Southern New Mexico. 
The main goal is to obtain photometry in five broad bands over 
10,000 deg? of high-latitude sky, and spectroscopy of a million 
galaxies and 100,000 quasars over this same region. SDSS 
DR7 comprises images and spectra for millions of galaxies, 
centered at a redshift ~0.1. The extinction-corrected spectro- 
scopic data are complete at 17.72 mag in the r-band. 

In previous studies, it was discovered that the SDSS 
photometric pipeline tends to overestimate the sky background 
of galaxies with extended low surface brightness by ~0.5 mag 
(Lisker et al. 2007; Liu et al. 2008; He et al. 2013). This is 
particularly relevant for dwarf galaxies with low surface 
brightness as their luminosities are more vulnerable to sky 
subtraction. To address this issue, we selected a sub-sample of 
low surface brightness galaxies (Du et al. 2015) in a.40, with 
central surface brightness in the B band j44(B)» 22.5 mag 
arcsec 7, and recalculated their photometry using the g-band 
and r-band SDSS mosaic images from the Twelfth Data 
Releases of the Sloan Digital Sky Survey (SDSS DR12; Alam 
et al. 2015), implementing improvements in sky-background 
subtraction following Zheng et al. (2015) and Guo et al. (2020). 
For galaxies with brighter surface brightness, we utilized 
Petrosian magnitudes from the released SDSS DR7 catalog. 

The stellar mass is estimated using the stellar mass-to-light 
ratio formula proposed by Bell et al. (2003) 


logig(M/L) = ax + by x (g — r). (1) 


Here we adopt the Kroupa stellar initial mass function 
(Kroupa 2002) with a,=—0.306—0.15 and b,= 1.097, 
a; = —0.222 — 0.15 and b; = 0.864 (Bell et al. 2003). 

The ALFALFA survey is the biggest blind extragalactic HI 
survey, with a wide sky coverage of 7000 deg’. The spatial 
resolution beam size is about 3/5. The catalog encompasses 21 cm 
line spectra of more than 30,000 extragalactic sources with radial 
velocities «18,000 km s !'. The average integration time of 
ALFALFA sources is 48 s, corresponding to a minimum flux of 
0.18 Jy at a velocity width of 50 kms~'. Here we use the 40% 
complete ALFALFA, a.k.a. o.40. There are 15,855 H I detected 
sources in a.40, among which 15,041 are extragalactic objects and 
the rest 814 are possibly to be Galactic high-velocity clouds. 
Cross-matched with the SDSS DR7, there are 12,423 galaxies 
with optical counterparts (Haynes et al. 2011). 

We select dwarf galaxies following the criteria in Hu et al. 
(2023). First, we select galaxies with absolute magnitude 
M,» —18 and the HI spectra with high signal-to-noise ratios, 
S/Ru;» 10. Then, we visually remove the galaxies which 
have multiple optical counterparts within 6’ (corresponding to 
twice the Arecibo beam size), which leaves 770 dwarf galaxies. 
Second, we use the regions where both SDSS and ALFALFA 
have detection, as shown in Figure 1. In order to avoid the 
incompleteness of neighbor galaxies at the edges of the survey, 
we further restrict to 5? < decl. < 35°, 125? < R.A. « 235°: 
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Figure 1. The footprint of surveys. The gray and green dots represent the SDSS spectroscopic galaxies within the redshift range of 0 < z « 0.017 and ALFALFA 
galaxies within the redshift range of 0.003 < z « 0.0133, respectively. The red rectangular regions represent the regions selected for the parent a.40 dwarf galaxies. 


Further details can be found in the main text. 


and 0° « decl. < 30°, —25? « R.A. « 25°, and we request the 
samples to be within a redshift range (0.003, 0.0133), leading 
to 413 selected galaxies. Finally, we restrict our sample to have 
r-band magnitude brighter than 20 and HI mass 105 ^M... 
This r-band limit is selected to have an optically complete 
sample set at 9596 level, and this HI mass limit is chosen to 
ensure an HI mass complete sample in the redshift of interest 
(Figure 2). The final sample consists of 320 dwarf galaxies. 


2.2. Deep Observations with FAST 


The ALFALFA survey has provided a substantial sample of 
galaxies with detectable HI content. However, due to the 
relatively short integration time of the survey, the sample 
mainly consists of H I-rich dwarf galaxies. In order to gain a 
comprehensive understanding of the variability of H I content 
in different environments, including H I-poor galaxies, we have 
increased the integration. time and utilized the FAST 
(Nan 2006; Jiang et al. 2019) to observe a set of 10 randomly 
selected dwarf galaxies, which have angular distances to any 
a.100 sources greater than 6’ (twice the Arecibo beam size). 
These galaxies either have not been detected in the ALFALFA 
survey or have low signal-to-noise ratios (S/Ns, as shown in 
Figure 2). By expanding our observations to include these 
targets, we aim to obtain a more comprehensive understanding 
of the distribution of H I content across various environments. 

FAST is a single-dish radio telescope with an effective 
diameter of 300 m and the highest sensitivity among ground- 
based single-dish radio telescopes. FAST is currently mounted 
with a 19-beam receiver covering a frequency range of 
1.05-1.45 GHz. Its frequency resolution is about 7.63 kHz 
(65,536 channels over 1.0-1.5 GHz bandwidth) and the angular 
resolution is about 2/9 at low redshift for 21 cm line spectra. 
The system temperature is ~20 K for zenith angles within 
26.4 deg observations. 


5 From COMBO S11 field, see on SDSS website: https: / /www.sdss4.org/ 
dr17 /imaging/other info. 


The observations were carried out based on a "Shared-Risk" 
project, using the position switch ON-OFF mode recording two 
polarizations (XX and YY) for each target. The ON-source and 
OFF-source integration time is 120s in most cases, and the 
overhead time (switch between ON-source and OFF-source 
positions) is 30s. We have carefully chosen the OFF point for 
each source so that a side beam (M08 or M14’) would be pointing 
to the source when the central beam (M01) is pointing to the OFF 
point (Zheng et al. 2020). This strategy successfully reduced the 
overhead time rates down to around 20%. For most of our targets, 
the sampling time is 0.1 s, we inject the high (10 K) noise diode 
signal for 0.1s every 1s during the observation. For AGC 
125644, we use a 0.5 s sampling time and inject the high noise 
diode signal for 2s every 20s. The majority of the observations 
have a total on-source integration time of approximately 
19.2 minutes. Details are listed in Table 1. 

We perform the following procedures for each beam and 
their corresponding polarizations (XX and YY). 

Flux calibration. We use the high (10 K) noise diode other 
than the standard calibrators for our flux calibrations. First, we 
convert the recorded noise diode signal into antenna temper- 
ature (K). To obtain the flux density in units of Janskys (Jy), we 
utilized the gain values reported in Table 5 of Jiang et al. 
(2020), for the M01 and M08 (M14) beams, respectively. 

Baseline and standing wave removal. Following Y. Jing et al. 
(2023, in preparation)? we remove a lower-order Polynomial 
baseline to each OFF-subtracted spectrum (ON-source subtract 
OFF-source for each beam, i.e., M01 ON-source position minus 
M01 OFF-source position?). There are ~1 MHz standing waves 


7 We found that some delta RFI only appeared in the M14 beam so we only 


used the M14 beam for those MOS positions with contamination. 
* HiFAST: https: / [ hifast.readthedocs.io. 


? Initially, our plan was to subtract the continuum using the same time from 
OFF position, specifically the M01 ON-source position minus the M08 ON- 
source position (pointing the OFF-source position). However, we discovered 
that the period and phase of the standing wave differed between the different 
beams, making it challenging to fit the standing wave if we performed the same 
time OFF subtraction. We noticed that the standing wave varies with time, but 
we were able to subtract the ripple of the spectra by performing the same beam 
OFF subtraction. 
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Table 1 
Observational Information 
Name R.A. Decl. z Distance Date tint/ trot Beam Comment 
degree(°) degree(°) (Mpc) (minutes /minutes) 
IC 3331 186.52215 11.81221 0.00411 17.6 6/30/2019 16/24 1, 14 
VCC 488 185.41364 7.25363 0.00921 39.5 6/30/2019 19.2/24 1,8 
IC 3487 188.30599 9.39735 0.00359 16.2 6/30/2019 19.2/24 1,8 No signal 
IC 3663 190.41421 12.24749 0.00309 15.7 7/18/2019 19.2/24 1,8 No signal 
UGC 7854 190.46653 9.40295 0.00324 15.9 7/18/2019 19.2/24 1,8 No signal 
LEDA 1718355 175.34892 24.85794 0.01280 54.9 8/19/2019 19.2/24 1,8 
LEDA 4477763 236.80355 17.95755 0.01110 47.6 8/19/2019 19.2/24 1, 14 
ECO 915 237.41397 18.88038 0.01328 56.9 8/19/2019 19.2/24 1,8 
UGC 8838 208.91125 4.98472 0.00390 22.2 8/19/2019 8/24 1 
AGC 125644 33.305 4.63583 0.02278 93.8 1/11/2022 15.7/20 1,8 PID: ZD2021_4 


Note. Columns: (1) Galaxy name, (2) and (3) coordinates, (4) spectroscopic redshift, (5) distance, (6) observed date, (7) integration time over total observing time unit 
in minutes, (8) FAST beam numbers which were used for integration, (9) comments. Most of our targets are observed during 2019 (PID: 2019a-003S), while AGC 


125644 is observed using the PID: ZD2021_4. 
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Figure 2. HI mass vs. redshift. Gray dots represent galaxies from the a.100 
sample, while the blue dots depict our volume-limited dwarf sample selected 
from a.40. Red stars correspond to the galaxies observed by FAST. The arrows 
denote the three FAST targets which have an upper limit of My ı < 107 M... 


Red dashed lines denote the selection boundaries. The black curve illustrates 


the detection limits of HI mass, assuming a velocity width of 50 km s~’. 


covering all bandwidth of FAST spectra, which are caused by 
reflections between the receiver and the bottom panel of the 
telescope (138.63 m causing a 1.08 MHz standing wave). We 
remove the standing wave by fitting a sinusoidal function using 
spectra around the heliocentric H I velocity of the target. 

The observations were conducted in topocentric mode, 
which means that the reference frame is not a rest frame due to 


the Earth's movement. To account for this, we convert the 
observed frequency to heliocentric velocity. The final spectra 
are obtained by averaging the spectra from MOI and 
M08 (M14). 

The spectra of detected FAST sources are presented in 
Figure 3. Out of these sources, seven of them have good S/Ns, 
which allows for reliable analysis of their HI properties. 
However, for the remaining three sources that have low S/N, 
their spectra only provide upper limits for the analysis. These 
upper limits are used to estimate the potential H I content in 
these sources and we incorporate them into the overall analysis. 

To illustrate the advantage of using longer integration time 
with the FAST telescope for HI, we choose two galaxies for 
comparison from the ALFALFA sample: UGC 8838 and AGC 
125644. The spectral profile of UGC 8838 in Figure 4 exhibits 
a relatively low S/N, while AGC 125644 displays a dubious 
sub-component between 6900 and 7000 kms !, also with a 
very low S/N. Figure 4 illustrates that, by utilizing a longer 
integration time, FAST effectively improved the S/N of UGC 
8838 from 7.5 to 14.8. For AGC 125644, we confirmed the 
existence of the suspicious sub-component, despite the S/N in 
ALFALFA already being as high as 7.1. This discovery led to 
an expansion of its line width by a factor of two compared to its 
measurement in ALFALFA. Overall, it demonstrates the 
improved detection capabilities of the FAST telescope through 
longer integration times, enabling the identification of pre- 
viously unidentified sub-components and enhancing the S/N 
ratio of the observed galaxies. 

We future extract the physical properties as follows, with the 
derived properties listed in Table 2. 

Velocity width Wsy and Wọ. The algorithm follows the 
procedure in Guo et al. (2020). We distinguish between 
systems exhibiting one peak and two peaks in their HI spectra. 
For sources with a single peak, we identify the wavelengths at 
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Figure 3. H I spectra of FAST targets. Each panel shows the spectrum of an individual galaxy detected by FAST using a long integration time (Table 1). Red curves 
are the FAST observations with a frequency resolution of 7.63 kHz (~1.65 km s~'), and blue curves depict the spectra smoothed using a Hanning smoothing window. 
The yellow regions indicate the 1o,ms regions of baselines. Blue dashed lines represent their heliocentric velocity (center velocity of Wso). 
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Table 2 
FAST Sources Parameters 

Name Wso Wao Vnelio log M, log Mui Rms S/N S/N Env Env 

(km s~') (km s~') (km s~') (Ms) (Mo) (mJy) Average ALFALFA Group N 
IC 3331 29.2 49.93 1237.6 8.65 7.09 0.73 7.9 9.6 1.03 56 
VCC 488 34.43 45.52 2772.0 8.35 7.19 0.59 8.3 10.9 7.56 19 
IC 3487 1074.8 8.34 6.40 | 0.66 0.93 68 
IC 3663 927.1 8.55 6.32 | 0.63 1.76 43 
UGC 7854 971.9 8.79 6.26 | 0.51 1.57 41 
LEDA 1 718 355 56.19 60.75 3861.9 8.57 8.23 0.62 6.9 11.5 3.16 6 
LEDA 4 477 763 79.07 106.72 3356.6 8.16 8.66 0.61 17.6 35 7.95 5 
ECO 915 67.83 89.41 4011.7 8.42 9.06 0.61 35.9 66.1 1.37 6 
UGC 8838 73.07 98.99 1182.7 9.04 7.79 0.94 TT 14.8 0.64 7 
AGC 125644 139.97 166.79 6872.1 8.88 9.12 0.55 8.2 21.8 8.56 


Note. Columns: (1) Galaxy name, (2) and (3) velocity widths, (4) heliocentric velocity, (5) stellar mass (from r-band estimation) (6) neutral hydrogen mass, (7) rms 
under 1.65 km '' velocity resolution, (8) and (9) signal-to-noise ratios, (10) distance to the nearest group environment, (11) projected number density environment. 
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Figure 4. Comparison of H I spectra between deep FAST observation and 
a.40. The blue dashed lines correspond to their heliocentric velocity (center 
velocity of Wso). The black curves depict the spectra from ALFALFA, while 
the red curves are the FAST observations with a frequency resolution of 
7.63 kHz (71.65 km s !). The solid blue curves depict the FAST spectra 
smoothed using a Hanning smoothing window. The yellow regions indicate the 
log, regions of baselines. The upper panel displays the H I spectrum of UGC 
8838, while the lower panel shows the spectrum of AGC 125644. 


which the fluxes are 50% (20%) of the peak values to obtain the 
raw values, W50,raw (W20,raw). For sources with two peaks, the 
50% (20%) fluxes are calculated for each side, and we identify 
the corresponding characteristic wavelengths at each rising side 
of the spectrum to measure the velocity width. We perform 
polynomial fitting at each rising side (Springob et al. 2005; 
Haynes et al. 2011) to obtain the characteristic wavelengths. 


We then subtract the instrumental broadening as follows 
(Catinella et al. 2012): 


Wsooo) = 4WA, — Av? (2) 


where Av represents the velocity resolution, which is 
approximately 1.65 km s 

The heliocentric velocity and H I line integrated flux density, 
S21, in units of Jy km s~'. The heliocentric velocity is defined 
as the center velocity of Wso. The total flux is then computed by 
integrating the continuum subtracted HI spectrum within 
the +2 x Ws 9 range from the heliocentric velocity. The 5c 
upper limits of the flux density for three undetected 
galaxies (IC 3487, IC 3663, and UGC 7854) are calculated 
as 54/Wé9/1.65 x o2,. x 1.65 under the assumption that 
Wso = 100 km s™'. 

Signal-to-noise ratio (S/N) of the spectrum. We perform two 
different estimations of S/N: 


(1) An S/N estimation analog to ALFALFA following 
Haynes et al. (2011), the S/N of the detection is 
estimated as: 


S/N = (S21/ W50) x whl? /onms (3) 


where wsmo is a smoothing width set to be Wso/20 for 
Wso < 400 km s™', and to be 20 for Wsy > 400 km s^, 
and Cms is the root-mean-square of the baseline. 

(2) An average S/N within velocity width Wso: 


(S/N)avg = (S21/Ws0) /Orms- (4) 


H mass and distance. In ALFALFA, the distances for galaxies 
with cz > 6000 km s~! are the Hubble distances, estimated 
as d = cz/ Ho where Hy = 70 km s7! Mpc! is the Hubble 
constant, while for galaxies with cz< 6000 km s^, the 
distances are obtained with a combination of primary 
distances from the literature and secondary from the 
Tully—Fisher relation (TFR). For our newly observed FAST 
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Figure 5. Absolute r-band magnitude distribution. The blue-filled histogram 
represents the parent a.40 dwarf galaxies, while the black and cyan histograms 
correspond to the TNGSO dwarf and TNG50 gas-rich dwarf galaxies. 


targets, we only use the primary distances in the literature, 
as they are more accurate. When no accurate distance 
measurements are available in the literature, we use Hubble 
distance instead which is calculated from the heliocentric 
velocity of the HI spectrum. The HI mass of the galaxy is 
then derived from the integrated HI flux and distance 
following Roberts (1962), under the assumption that the H I 
content is optically thin: 


My; = 2.36 x 105d? f S21dv [Mo] (5) 


where d is the luminosity distance to the galaxy unit in Mpc. 


2.3. Dwarf Galaxies in TNG50 


We make use of the Illustris-TNG cosmological hydro- 
dynamical simulations? for the comparison. TNG50 (Nelson 
et al. 2019; Pillepich et al. 2019) traces 2 x 2160? dark matter 
particles and gas cells in a period box size of 51.7 Mpc. The 
mass resolution is 4.5 x 10°M, and 8.5 x 10‘M, for dark 
matter and baryon, respectively. Stellar masses are calculated 
using all star particles within each subhalo. We adopt the 
"Neutral Hydrogen Abundance" from the catalog for H I mass 
in non-star-forming cells and calculate the cold hydrogen gas 
mass (HI J-H5) using the two-phase interstellar medium model 
(a modification of Springel & Hernquist 2003) in star-forming 
cells. Then, we drop out the molecular hydrogen from all cold 


1) https:/ /www.tng-project.org 
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Figure 6. Environment distributions. The left panel shows the histogram of the 
distance to the nearest group, while the right panel displays the histogram of the 
projected number density. The blue histograms represent the parent a.40 dwarf 
galaxies, while the red histograms correspond to the FAST targets. The black 
vertical solid lines indicate three times the virial radius (left panel) and the 
median number density value (right panel) of the a.40 dwarf galaxies. 


hydrogen mass to obtain the H I mass for star-forming cells. 
The molecular hydrogen fraction fy, in star-forming regions is 
calculated by the modified KMT model (Krumholz et al. 2009; 
McKee & Krumholz 2010) as 


_fl-3s/4+s) ifs<2 
fa = E ifs >2 (6) 


where s is given by 


sa Ind_+.0.6x + 001) 


7 
0.67, 7) 
and 
x = 0.756(1 + 3.179365) (8) 
LgasT 
Te = — (9) 
" 


where Z represents the metallicity for gas cell in units of solar 
metallicity Zo = 0.0127, o =Z x 10” m? is the dust cross- 
section and u = 3.9 x 10?" kg is the mean mass per particles. The 
Jeans length approximation method is adopted when estimating the 
gas surface density Xas. We discard galaxies with less than 200 
star particles to ensure a reliable estimation of stellar mass. 

In order to compare to observations, here we select dwarf 
galaxies with M, > —18 and HI mass 105 ^M. for the main 
comparison (TNG50 gas-rich dwarf galaxies). While we did 
not set a lower limit for the r-band luminosity, Figure 5 
demonstrates that the distributions in M, are similar between 
the a.40 dwarf sample and the TNG50 gas-rich dwarf galaxy 
sample. Furthermore, to examine the impact of environmental 
effects, we utilized a larger sample by including all galaxies 
with M, > —18, which encompasses both gas-rich and gas-poor 
galaxies. As depicted in Figure 5, it is evident that such 
selection also encompasses a broader range of faint galaxies in 
terms of r-band magnitude. 
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Figure 7. Left: H I mass vs. stellar mass relation. The red dashed line represents the completeness limit of My ; = 105! ^M... Right: the H I-to-stellar mass ratio vs. 
stellar mass relation. The black and cyan contours correspond to the TNG50 dwarf and TNG50 gas-rich dwarf samples, respectively. The blue dots represent parent 
a.40 dwarf galaxies, while the red stars indicate the FAST targets. The black lines depict the upper limit of the contours for the TNG50 sample. 


2.4. Environment 


We adopt two kinds of environment definitions, the projected 
number density and the distance to the nearest group/cluster. 

The projected number density is defined as the number of 
galaxies brighter than —16.6 mag (r-band) within 1 Mpc 
projected distance and up to +1000 kms ' velocity difference 
along the line-of-sight. Here we use the spectroscopic data of 
SDSS which is complete at 17.72 mag in the r-band, 
corresponding to an absolute magnitude of ~—16.6 mag at 


the redshift of 0.017 (z = 0.0133 + ks"), 

Distance to the nearest group is defined as the minimum value 
of the distance to the group center in a unit of the corresponding 
group viral radius (Guo et al. 2020). The group catalog is taken 
from (Saulder et al. 2016) which applies a friends-of-friends (FoF) 
group finder algorithm on the SDSS DR12 (Alam et al. 2015) and 
the 2MRS (Huchra et al. 2012). Compared to the group catalog by 
Yang et al. (2007), the group catalog used in this paper includes 
more groups in the local Universe (including groups at z « 0.01). 
Various observational biases have been taken into account, 
including the Malmquist bias, the “Fingers of God," etc. The 
groups/clusters masses are derived from the total luminosity, 
luminosity distance, velocity dispersion, group radius, and the 
number of detected group members (we refer the reader to Saulder 
et al. 2016 for more details). This method has been precisely 
calibrated on the mock catalogs from the Millennium simulation 


(Guo et al. 2011, 2013). The minimum and median distances to the 
nearest group of the parent dwarf galaxies are 0.074 and 4.339, 
respectively. 

The histograms of two different definitions of the environment 
are shown in Figure 6. Most a.40 dwarf galaxies reside in 
relatively low-density regions, outside the virial radius of 
surrounding clusters/groups. Galaxies at high densities or close 
to groups/clusters suffer from environmental stripping and likely 
lose their HI gas in a short timescale. One needs a longer 
integration time to reach deep detection of galaxies in such 
environments. FAST, with its longer integration time, captures a 
higher fraction of galaxies in these high-density regions, as 
depicted by the red histogram. Accordingly, the combination of 
the gas-rich a.40 galaxies and the dwarf galaxies residing in 
denser regions, exposed for a longer duration through FAST, form 
a valuable sample for investigating environmental dependencies. 


3. Results 


In this section, we investigate the relationship between the 
H I and stellar mass, their environmental dependence, and the 
scaling relations between the HlI-to-stellar mass ratio and 
optical properties. 


3.1. HI to Stellar Mass Relation 


Figure 7 presents the relationship between HI mass and 
stellar mass. The H I mass of the jointly selected gas-rich dwarf 
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Figure 8. H I-to-stellar mass ratio as a function of environments. Upper left panel: Results for the observational data. The blue dots are the a.40 dwarf galaxies, with 
the red line denoting the median My ;/M,, value. The orange and the green dots are dwarf galaxies from SHIELD and LITTLE THINGS. The red stars are our FAST 
targets, and the arrows indicate the upper limit of undetected sources. The black line indicates the critical radius of three times the virial radius. The environment is 
defined by the distance to the nearest group. Upper right panel: Similar to the upper left panel but the environment is defined by projected number density. The black 
line represents the median number density. Lower left panel: similar to the upper left panel but for the TNG50 galaxies. The gray and cyan dots represent the TNG50 
dwarf and TNG50 gas-rich dwarf galaxies. The dark and blue curves represent the corresponding 16%, 50%, and 84% percentage levels. Lower right panel: Similar to 
the upper right panels but for the TNG50 galaxies. Data have the same line style and color coding as those in the lower left panel. 


galaxy sample appears to increase with their stellar masses. The 
observed H I versus stellar mass relation is well reproduced in 
simulations when taking into account the detection limit of the 
a.40. However, the gas-poor galaxies detected by FAST push 


this H I mass versus stellar mass relation downwards more than 
an order. 


When including both gas-rich and gas-poor dwarf galaxies in 
simulations, it shows an even stronger dependence of the HI 
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mass with stellar masses. Dwarf galaxies with low H I masses 
as detected by FAST are also found in simulations. Indeed, 
TNG50 exhibits a continuous distribution along the y-axis, with 
a relatively small fraction of galaxies having log My, < 7.5 
at log M. ~ 8.5. 

Interestingly, we observe a common boundary for the 
maximum Mg ; (HI mass) in both our observational data and 
simulations. This boundary sets a limit on the highest H I mass 
that can exist for a given stellar mass (M). Such behavior also 
translates into the My; to M. ratio as a function of stellar mass 
(right panel). To describe this upper boundary, we employed a 
fitting function as follows 


logig(Mui/M&) = —0.5 x logig(Mx) + 5.2. (10) 
This is equivalent to 
2 
Min 10104... (11) 
* 


Given the detection limit, at least 97.2% 0.40 dwarf galaxies 
are below this boundary, while in simulations for TNG50 dwarf 
galaxies, 99.296 are below this boundary. 


3.2. H r-to-stellar Mass Ratio in Different Environments 


Figure 8 displays the distribution of the H I-to-stellar mass 
ratio (My j/ M.) across different environments. To extend the 
sample size, we include data from two Local Volume dwarf 
galaxy samples: Local Irregulars That Trace Luminosity 
Extremes, The HI Nearby Galaxy Survey (LITTLE THINGS; 
Oh et al. 2015) and Survey of HI in Extremely Low-mass 
Dwarfs (SHIELD; McNichols et al. 2016). 

The left panels of the figure demonstrate a slight decrease in the 
My 1/ M,, as the distance to the nearest group decreases. The black 
lines in the left panels indicate where the distance to the nearest 
group is three times the group virial radius, beyond which the 
high-density environmental effects typically vanish. In the right 
panels, the black lines correspond to the median values of the 
number density for the a.40 dwarf galaxies (3) and TNG50 dwarf 
(1). A consistent trend is observed, indicating a decrease in the 
My 1/ M. as the environment becomes denser, regardless of the 
specific definition of the environment. To remove the mass 
dependence, Figure 9 shows the My ;/M,, versus M. relations for 
galaxies in the 25% most dense and the 25% most isolated 
environments. It is clear that galaxies in dense regions have a 
lower My 1/M,, at any given stellar mass. 

We observe a lower boundary of the My ;/M,, ratio at around 
0.1, which could be an artifact resulting from the limited 
observation depth of the a.40 survey. At low densities galaxies 
detected by FAST reside in the same region as the a.40 sample. 
However, the high sensitivity of FAST, along with longer 
integration times, extends the lower boundary by a factor of 10 
near groups/clusters. These findings suggest that there is a 


10 


Hu et al. 


dist to nearest group number density 


log Muj/M« 
log Muj//M« 


* ad0 field 
*  a40 group 


* ad0 field 
*  a40 group 


-2 2 
65 70 75 80 85 90 65 70 75 


log M4[MG] 


95 100 105 85 90 95 100 105 


80 
log M4[MG] 


Figure 9. H I-to-stellar mass ratio as a function of stellar mass in different 
environments. Left panel: The environment is defined by the distance to the 
nearest group. Right panel: The environment is defined by the number density. 
Blue dots are the a.40 dwarf galaxies in the 25% most isolated environments, 
while red dots are galaxies in the 25% most dense environments. 


gradual and continuous process of HI gas stripping around 
groups/clusters, rather than a sudden removal. 

Such environmental dependence is much weaker, if any, in 
simulations. In the case of TNG50 gas-rich dwarf galaxies 
which mimic the observed a.40 dwarf sample, their My 1/Mx 
ratio does not vary much when galaxies get closer to groups/ 
clusters, or when the density increases. This suggests that gas 
stripping occurs faster in simulations than those in the real 
Universe. When including gas-poor dwarf galaxies in TNG50, 
it encloses the regions covered by the gas-poor dwarf galaxies 
detected through the deep FAST observation. The dependence 
on environments is stronger compared to the TNG50 gas-rich 
sample, but overall, the dependence remains weak. 

In summary, the process of gas stripping from galaxies as 
they approach high-density regions is rapid, but the timescale 
of this process could be even shorter in simulations than those 
in observations. 


3.3. H I-to-stellar Mass Ratio versus Galaxy Optical 
Properties 


Previous studies found that H I-to-stellar mass ratio is closely 
related to galaxies' optical properties (Zhang et al. 2009; Zhang 
et al. 2021; Catinella et al. 2010; Li et al. 2012). In Figure 10, we 
present the relation between My ;/M,, and optical proprieties, as 
proposed by Zhang et al. (2009). In the upper left panel, we show 
the My 1/ M, as a function of the combination of the g — r color 
and stellar surface mass density. The stellar surface mass density is 
defined as log(u, ;) = log(M&) — log(27R%,;), where M, 
represents the total stellar mass, and R50; corresponds to the i-band 
Petrosian half-light radius. In the upper right panel, the stellar 
surface mass density is replaced with surface brightness, which is 
defined as u; = m; + 2.5 logQzRs, ). where m; represents the 
apparent Petrosian i-band magnitude. Zhang et al. (2009) used a 
sample of 800 galaxies from the SDSS DR4 and the HyperLeda 
H I (Paturel et al. 2003) catalog (Sample II) and obtain the scaling 
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Figure 10. Scaling relations between H I to stellar mass ratio and optical properties. Upper left panel: My 1/ M. vs. f(g — r, Hx) . The dots and stars represent the a.40 
dwarf galaxies and the FAST targets. Symbols are colored by log(number density +1). Upper right panel: Similar to the upper left panel but for My 1/ M, vs. f(g — r, 
Hi). Lower left panel: Similar to the upper left panel but for TNGSO gas-rich dwarf galaxies. Lower right panel: Similar to the upper left panel but for TNG50 dwarf 
galaxies. The black solid lines and dashed lines indicate the scaling relation and its 1c scatter established by Zhang et al. (2009). 


relations as presented by the solid lines, with the dashed lines 
corresponding to the lo scatter. 

Figure 10 demonstrates that the H I-selected sample from 
a.40 exhibits a slightly higher Mu1/M, compared to the 
sample studied by Zhang et al. (2009). This divergence could 
be attributed to the exclusion of gas-poor galaxies in a.40, 
which are primarily found in high-density regions. To address 


11 


this, we incorporate the inclusion of high-density gas-poor 
galaxies detected by the FAST telescope, particularly in high- 
density environments. When analyzing the samples separately 
based on different environments, a clear dependence on the 
environment becomes apparent. For a given combination of the 


optical properties, galaxies in dense regions seem to have lower 
My 1/My. 
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In the lower panels of Figure 10, we present a comparison 
with the TNG50 gas-rich dwarf (left) and TNG50 dwarf 
galaxies (right). For gas-rich dwarf galaxies, it shows a 
dependence on environments mainly on the right-hand side 
where the gas fraction is relatively high. When including gas- 
poor dwarf galaxies, My ;/M,, reaches lower values at the left- 
hand side, galaxies in dense regions are presented well below 
the scaling relations found by Zhang et al. (2009). 

These findings imply that the current understanding of the 
relationship between the My ,/M,, and optical properties may 
not be accurate enough without considering the influence of 
environmental factors. 


4. Summary 


In this paper, we investigate the H I-to-stellar mass ratio 
(My 1/ M), the scaling relations between My ;/M,, and optical 
properties, and their dependence on environments using a 
jointly selected volume-limited sample from a.40 survey and 
SDSS. Our sample primarily consists of gas-rich dwarf 
galaxies, and to complement this data set with more gas-poor 
galaxies, we incorporate data from 10 dwarf galaxies observed 
with the FAST telescope. These additional observations were 
conducted with significantly longer integration times, allowing 
for a more comprehensive analysis. Our main findings are as 
follows: 


(1) There exists a clear correlation between HI mass and 
stellar mass in our jointly selected sample of dwarf 
galaxies. This finding is also supported by simulations. 
Furthermore, we observe a common upper boundary for 
HI mass among both the simulation and observational 
data sets for each given stellar mass. This suggests that 
there is a limit to the amount of HI mass that dwarf 
galaxies can possess, regardless of their stellar mass. 

(2) The H I-to-stellar mass ratio exhibits a decreasing trend as 
the environment becomes denser, regardless of whether it 
is defined by the distance to the nearest group or the 
(projected) number density. These findings are further 
supported by our analysis of deeper observations of high- 
density gas-poor targets. Simulations also exhibit similar 
trends, but the slope is much shallower. This suggests that 
the stripping of gas in high-density regions occurs more 
gradually in the real Universe compared to simulations. 

(3) The a.40 dwarf galaxies show a slight departure from the 
median of the My ;/ M, versus optical properties scaling 
relation established by Zhang et al. (2009), with a 
tendency toward higher Myj;/M, values while still 
mostly falling within the one-sigma region. We also 
observe a dependency on the environment for both the 
observed and simulated data sets, although this effect is 
less pronounced in the simulations. 
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The majority of our findings are derived from the a.40 dwarf 
galaxies. Although we have added a gas-poor sample with 10 
dwarf galaxies through deep observations with FAST, we 
acknowledge that the statistical significance remains limited. 
Therefore, a more comprehensive and extensive survey is 
necessary to determine if these galaxies are indeed outliers or if 
they simply contribute to the scatter observed in the scaling 
relations. 
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